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The finite temperature lattice QCD with Nf = 2 nonperturbatively improved Wilson 
*~ j" 1 ■ fermions is studied on 16 3 8 lattice. Using abelian projection after fixing to MA gauge we 
^-i. determine the transition temperature for m n /m p ~ 0.8. 
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1. Introduction 

- T— I | 

^ ■ Determination of the critical temperature of the chiral phase transition is one of the im- 
portant nonperturbative problems in QCD to be addressed in lattice simulations. Bielefeld 
group and CP-PACS collaborations using two different types of improved lattice action 
for fermions were able to estimate T c in the chiral limit and their values are in good agree- 
ment j 1 112 j . Still there are many sources of systematic uncertainties and new computations 
of T c with different actions are useful as an additional check. We made first large scale 
simulations of the nonperturbatively 0(a) improved Wilson fermion action. Moreover we 
performed simulations with the lattice spacing a substantially smaller than in studies by 
Bielefeld and CP-PACS groups. 

The following fermionic action is employed in our study: 

S F = sP - ]-Kgc sw a b ^ti)(x)a^ v F^(x), (1) 

^ X 

where Sp is the original Wilson action, c sw is the clover coefficient determined nonper- 
turbatively jS] • We use Wilson gauge field action. 

The action has been used in T = studies of lattice QCD by UKQCD and QCDSF 
collaborations |4|5j . These studies confirmed that 0(a) lattice artifacts are suppressed as 
expected. To fix the physical scale and m^/rrip ratio we used results obtained by these 
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Table. 1 Simulation statistics. 



collaborations 0J. So far only N t = 4 and 6 finite temperature results obtained with 
action (0) are available |JJ. These results were obtained at a rather large quark mass 
(m^/rrip > 0.85). In this work me make simulations with N t = 8 what allows us to 
decrease m n /m p down to ~ 0.8. We choose the spatial extension of the lattice N s = 16 
as a compromise between computational burden and need to reduce the finite size effects. 

It is known that the order parameter of the finite temperature phase transition in 
quenched QCD is the Polyakov loop and the corresponding symmetry is global Z(3) 
symmetry. In the chiral QCD the order parameter of the chiral symmetry breaking 
transition is the chiral condensate < ipip >. As numerical results show PQ both order 
parameters can be used to locate the transition point at intermediate values of the quark 
mass. We use Polyakov loop susceptibility to determine the transition temperature. 

2. Simulation details 

Hybrid Monte Carlo algorithm with parameters 5t = 0.0125, n T = 20, providing accep- 
tance rates of about 70% is used in our simulations. The simulations were done on SR8000 
Hitachi at KEK, Tsukuba and MVS 1000M at Joint Supercomputer Center, Moscow. For 
analysis SX5 NEC at RCNP and PC-cluster at ITP, Kanazawa were employed. Our code 
performs at the speed of 2.4 GFlops per node on Hitachi computer. We needed from 1000 
(r = N tra j -n T ■ 5t = 250) to 3000 (r = 750) trajectories for thermalization, depending on 
k and f3. For runs started from configurations generated at the adjacent k this value was 
much lower. We determined the transition temperature at two values of (3, 5.2 and 5.25, 
varying k as shown in Table 1. 

We fixed the maximally abelian (MA) gauge on generated configurations. The simulated 
annealing algorithm has been used to fix the gauge. The advantages of this algorithm 
in comparison with the usual iterative algorithm has been demonstrated in SU(2) pure 
gauge theory for MA gauge |S] and Maximal Center gauge p. 

Abelian and monopole Polyakov loops and their susceptibilities were measured on gauge 
fixed configurations. We found that for abelian and monopole observables the signal/noise 
ratio was better than that for gauge invariant nonabelian observables. This observation 
is in agreement with the results from quenched QCD at T > and both quenched and 
un quenched QCD at T = 0. In particular the maximum of the Polyakov loop susceptibility 
is essentially better separated from the rest of the data for the monopole Polyakov loop 
than for nonabelian Polyakov loop. It is important that the transition temperature values 
determined by both susceptibilities are the same. 
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Figure 1. Nonabelian Polyakov loop at both /3's (left), and the abelian, monopole and photon Polyakov 
loops at = 5.2 (right). 



3. Results and Conclusions 

In FigCJwe show results for average of various kinds of Polyakov loops. One can see that 
(P) is a smooth function of k and that abelian and monopole Polyakov loops behavior 
is qualitatively the same as behavior of the nonabelian Polyakov loop while the photon 
Polyakov loop is almost constant across the transition. 

Transition point k, denoted as K t , has been determined from the maximum of the 
Polyakov loop susceptibility, see Fig [21 As can be seen from Fig |2] the susceptibilities for 
abelian and monopole Polyakov loops have maxima at the same value of k. We found 
that Kt = 0.1344(1) at (3 = 5.2 and K t = 0.1341(1) at (5 = 5.25, see FigH This was 
transformed into transition temperature with the help of the interpolation formula for 
r /a [4 . The results are: T c r = 0.54(2) and 0.56(2). In physical units, taking l/r = 394 
MeV, we obtained T c = 213(10) and 222(10)MeV, respectively. Using again data from [1] 
we estimate m^/rrip = 0.78, 0.82 at the transition points. 

In FigElwe show our results for transition temperature in comparison with results of 
refs. |I] and 0. The ratio T c j \fo was computed for y/a = 425 MeV used in pp. We 
conclude that our results are in good qualitative agreement with previous results 0. 
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Figure 2. Nonabelian Polyakov loop susceptibilities at both f3's (left), abelian and monopole Polyakov 
loop susceptibilities at (3 — 5.2 (right) 
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Figure 3. Transition temperature 



This agreement implies that the dependence of the transition temperature on the lattice 
spacing is rather weak. 

For our lattice with N s /N t = 2 the question of the finite volume effects is very impor- 
tant. To check this effect we made simulations at (3 — 5.2, k = 0.1343 on iV s = 24 lattice. 
The k value was chosen close to the transition point where finite volume effects should 
be more pronounced. We found that both average of the Polyakov loop and its suscepti- 
bility are the same within error bars as in our main simulations on N s = 16 lattice. This 
implies that the finite volume effects should not spoil our conclusions. We are planning 
to perform simulations on 24 3 • 10 lattice at m n /m p < 0.7 to determine T c closer to the 
chiral limit. 
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